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Abstract 

We investigate the fluctuations around the average density profile 
in the weakly asymmetric exclusion process with open boundaries in 
the steady state. We show that these fluctuations are given, in the 
macroscopic limit, by a centered Gaussian fleld and we compute ex- 
plicitly its covariance function. We use two approaches. The flrst 
method is dynamical and based on fluctuations around the hydrody- 
namic limit. We prove that the density fluctuations evolve macro- 
scopically according to an autonomous stochastic equation, and we 
search for the stationary distribution of this evolution. The second 
approach, which is based on a representation of the steady state as 
a sum over paths, allows one to write the density fluctuations in the 
steady state as a sum over two independent processes, one of which 
is the derivative of a Brownian motion, the other one being related 
to a random path in a potential. 

Key words: Exclusion process, stationary non-equilibrium states, fluctu- 
ations 

We thank J. L. Lebowitz for useful discussions. 



*Laboratoire de Physique Statistique, Ecolc Normale Superieure, 24 rue Lhomond, 
75005 Paris, France; emails derrldaaips . ens . fr , enaudOlps . ens .f r 

tiMPA, Estrada Dona Castorina 110, CEP 22460 Rio de Janeiro, Brasil and 
CNRS UMR 6085, Universitc de Rouen, 76128 Mont Saint Aignan, France; e-mail 
landini@impa.br 

tCeremade, UMR CNRS 7534, Universite de Paris Dauphine, Place du Marcchal De 
Lattre De Tassigny 75775 Paris Cedex 16 - France; e-mail olla8ceremade.dauphine.fr 



1 



1 Introduction 



Non equilibrium systems such as systems in contact with two thermostats at 
unequal temperatures or with two reservoirs at unequal densities are known 
to exhibit long range correlations in their steady state 1 . These long range 
correlations have been calculated from the microscopic dynamics only in 
very few cases: mainly in the case of the symmetric exclusion process |21 
and of the asymmetric exclusion process [Sj . 

In the present paper we focus on the weakly asymmetric exclusion pro- 
cess (WASEP) (for which the bias scales as the inverse of the system size). 
We show that the fluctuations of density are Gaussian and that, as for the 
symmetric case, the direct calculation of the two point function by fluc- 
tuating hydrodynamics agrees with the expression derived by expanding 
the large deviation function around the average density profile. We also 
show that, as for the asymmetric case the correlation functions can be 
expressed in terms of two independent random processes. 

The asymmetric exclusion process describes the stochastic evolution of 
a system of particles on a one dimensional lattice of L sites. Each site of 
the lattice is occupied by at most one particle at a given time. Thus a 
configuration of the system is specified by a series of occupation numbers 
{7?i}j=i,- -,i: where 77^ = 1 if site i is occupied by a particle, 77^ = if site i 
is empty. 

At any given time, each particle independently attempts to jump to 
one of its neighboring sites on the lattice with rates which depend on the 
direction of the jump. The jump succeeds only if the target site is empty, 
otherwise the particle doesn't move. We choose to scale the time such that 
the hopping rate to a right neighbor is 1, and we note by q the hopping 
rate in the opposite direction. Physically, the asymmetry between left and 
right hopping rates mimics the effect of some external field acting on the 
particles. 

The first and last sites of the lattice, respectively labeled 1 and L, are 
in contact with reservoirs of density respectively pa and pb- This can be 
achieved by adding, to the left of site 1 , a site whose probability of being 
occupied by a particle is kept constant to pa independently of the rest of 
the system, and, to the right of site L, a site L -f 1 whose probability of 
being occupied by a particle is kept constant to p^ independently of the rest 
of the system. Another equivalent way of expressing the action of the left 
reservoir is to say that particles are added to the site 1 at rate a — pa when 
the first site is empty, and when it is occupied, the particle is removed from 
the system at 7 = q{l — pa) (i.e. a particle at site 1 attempts to jump to 
the left reservoir (site 0) with a rate q and the jump succeeds with a rate 
1 — Pa). In the same way, if site L is occupied, the particle is removed at 
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rate f3 — 1 — Pb, and if it is empty, a particle is added at rate 7 = qpb- 

The generator L of the dynamics, acting on a given function g of the 
configuration of the system 77 = {r?i}i=i,...,L is given by 



L-l 



9iv)] 



i=l 



L 



i=2 



+ [pa{i - m) + 9(1 - pa)m] W) - giv)] 
+ [Qpbii - til) + (1 - PbW] [giv"") - giv)] , 



where is the configuration obtained from i] by exchanging the occu- 

pation numbers of sites i and i + 1 and is obtained from 77 by changing 
the occupation number of site i. 

In the following, we consider the weakly asymmetric exclusion process, 
where the asymmetry q scales with the system size L by 



Of particular interest are the macroscopic properties of the system in 
the large L limit. We focus here on the macroscopic density profiles {p{x)}, 
< X < 1, obtained by rescaling by and smoothening the microscopic 
density profiles {77i}i=i,L. To do that, we may for example divide the 
system in mesoscopic boxes of size Lk, with 1 <^ Lk <^ L. The macroscopic 
density Pl{j^) is then simply the number of particles in the box containing 
the site i, divided by the size of the box. A more mathematical approach 
consists in defining, for each size L, the following distribution acting on 
smooth test function 



In spite of the stochastic evolution of the system at the microscopic scale, 
it is known that, as L — > 00, the macroscopic profile p^ converges almost 
surely to a deterministic evolution p^{t,x) given by the solution of the 
hydrodynamic equations of the process [HlEm. In the case of the weakly 
asymmetric exclusion process, it has been proved [S] [U] that this is given 
by the viscous Burgers equation: 




(1.1) 



dtp" = dip" - A(l - 2p")d^p 

p"{t,0)=pa, p"{t,l)^Pb, 

p"{0,x)^ po{x) 



(1.2) 
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where po{x) is the macroscopic Umit profile at time t — 0. 

At a finer scale, the density profile pL has random fluctuations of order 
L^^/^ around the hydrodynamic trajectory . One defines the density 
fluctuation field by 

iL{t.x)^^/L[pL{t,x)~ p"{t,x)\ . (1.3) 

It is known [Sj^j that these fields have a well defined large L limit C(x,t) 
which is a generalized Ornstein-Uhlenbeck process in the case of WASEP 
on an infinite lattice (for such system, the parameter L is not any more 
related to the size of the lattice, which is infinite, but only to the rescaling 
of time and space, as well as to the asymmetry rate). 

For a finite system of size L, after a long time, the system eventually 
reaches a steady state where the properties of the system become time- 
independent. Except for particular values of the reservoir densities and of 
the asymmetry parameter q for which detailed balance is satisfied JI] , this 
stationary state is a non-equilibrium steady state, which differs from an 
equilibrium state by the presence of a non-zero, site-independent average 
current j 

] = - Vi+i)) - - m)) (i-4) 

(where the brackets (.) stands for the average with respect to the steady 
state probability denoted by p^). 

The system presents almost surely an average macroscopic profile {p(x)} 
such that, for any site i, one has in the large L limit: 

Putting dtp to into (|1.2|l gives the following equation for the steady state 
average profile {p{x)}: 

p'{x)^Xp{x){l-p{x))-J, 

m = Pa , (1.5) 

= Pb , 

where J is an integration constant solution of 

j-P' dp 

xp{i -p)~j-^ • 

J is related to the steady state average current j (|1.4|l by 

J — lim Lj 

L — *oo 
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The steady state fluctuation field will be noted 



e\x)^ Mm VL[pLix) - pix)] 
with pl{x) defined as in ll.l|l . 

In the stationary state the macroscopic limit of the density profile coin- 
cides with the deterministic solution of the stationary equation (|1.5(l with 
the steady current J determined by the asymmetry parameter A and the 
boundary conditions Pa, pb- Actually, there exist ^Jj explicit expressions 
for but we will not use them in this paper. Without loss of generality 
we assume that 

Pa < Pb 

so that we have always p'{x) > 0. The current J can be positive or negative 
depending on A. 

Our main result is that the macroscopic fiuctuation field is a 

centered Gaussian field on [0, 1] with covariance given by 

= Xipi^mx -y) + f{x, y) (1.6) 

where 

X{p)=p(l-P) (1.7) 

and /(x, y) is given by 

/(„, . , (i,s, 

Jo p'{u) 

where x Ay = minja;, y} and xV y = maxjx, y}. 

In terms of the two points correlation function, this result implies that 
when X ^ y 

L {{v[Lx]VlLy]) - {V[LX]) {v[Ly])) f{x,y) (1.9) 

Some comments are in order. 

• In the symmetric case (A = 0) we have p'{x) = pt — Pa J = 
-{pb - Pa), so that 

f{x,y) = -{pi,- paf{xAy)[l-{xyy)] = ^{pt,^paf{-^)-\x,y) , 

in agreement with the result of Spohn on the 2-point correlation func- 
tion 12 . 
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• In the equilibrium case J = 0, the stationary fluctuation fleld is just 
white noise on [0,1], corresponding to the fact that the stationary 
state will be given by a product measure. 

• The sign of /(x, y) depends on the sign of the current J. If A is posi- 
tive and large enough, J > (the weak shock regime) and f{x, y) > 0. 
In the other cases J and / are negative. 

We derive the result H1.6I1.8|I by two different methods. The first ap- 
proach is dynamical: we search for the stationary solutions of the macro- 
scopic fluctuations process. The second approach is static and based on a 
representation (valid only when J{pa — Pb) > 0, see jTJj) of the weights in 
the steady state as sums over paths which was used in to calculate the 
large deviation function of the stationary measure. A priori there is no rea- 
son that the correlation functions are simply related to the large deviation 
functional of the density profile: density fluctuations describe variations 
of density of order whereas the large deviation functional describes 
variations of order 1. For the symmetric exclusion process, it was however 
shown |5j that the expression of correlation functions obtained by a direct 
calculation can be recovered by expanding the large deviation functional 
around the average profile. On the other hand, for the asymmetric ex- 
clusion process, the large deviation functional is non-analytic close to the 
average profile and there is no simple connection between the large devia- 
tions and the two point correlations T . One outcome of the present work 
is that the fluctuations of the density are still given, for the WASEP, as 
the expansion of the large deviation function around the average profile. 

Both methods rely on distinct representations of the kernel /. In the 
dynamical approach we proceed as follows. Let C be the differential oper- 
ator A + A(l — 2p(x))V with Dirichlet boundary conditions on [0, 1]. Since 
p is the solution of (|1.5|l . it is easy to show that C can be extended to a 
negative self-adjoint operator on {p' {x)dx) . The kernel of the inverse op- 
erator (— £)^^(a;,y) can be calculated explicitly (cf. Section|2l and shown 
to satisfy 

J{x,y) = J-p'{x)p'{y){-C)-\x,y) . (1.10) 

In the static approach, one shows that the Gaussian field £,^*{x) can be 
decomposed as a sum of two independent processes 

e{x)^sl^^B'{x) + \Y'{x), (1.11) 
where B{x) is a standard Brownian motion such that 

{[B{x)-B{x')f)^\x-x'\ (1.12) 
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(i.e. B'{x) is a standard 5-correlated white noise), while Y{x) is a centered 
Gaussian process whose distribution is formally given by 

,(1-13) 

where is the standard Feynman measure. Writing the distribution 

of the centered Gaussian process Y'ix) as 

dQ{{Y'}) cx exp |-i ^ Y'{u)T{u, v)Y'{v) du duj V[{Y'}] , (1.14) 

we show in section 4 that its covariance can be written as 

{Y'{x)Y'iy)) = T-\x, y) = 2x{p{x))5{x - y) + 4/(x, y) (1.15) 

where f{x,y) can be explicitly calculated and shown to be given by H1.8|l . 
The covariance (|1.6|l of the fluctuation field C'**(a;) follows from Hl.ll|l and 

2 Dynamical approach 

Consider the time dependent fluctuation field ^i,(i, x) defined by H1.3|l . 
Under proper assumptions on the initial distribution, it is proved in |T2| 
that the law of converges, as L — > oo, to the solution ^{t, x) of the linear 
stochastic partial differential equation: 

dt^ = dl(~XdA[l - 2p"{t,xm - (^^2x{p"{t,x)) W{t,x)^ , (2.1) 

where W{t, x) is the standard space time white noise, i.e., the Gaussian 
process on M.^ x (0,1) with covariance 

{W{t,q)W{t',q'))^Siq-q')Sit-t'), 

and X is given by H1.7|) . 

Since ^{t,x) and W{t,x) are distributions-valued processes, equation (|2.1() 
should be interpreted in the weak form: for any smooth test function G 
with compact support in (0, 1), denoting by ^t(G) = G{x)£,{t,x)dx, 

6(G)-eo(G) = f UG")ds + \ f t{[l-'ip"{s,-)]G')ds{2.2) 
Jo Jo 

+ f f G'{x)J2x{p"{s,x)) W{s,x)dsdx . 
Jo Jo 
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To investigate the asymptotic behavior of ^t, consider the weakly asym- 
metric exclusion process starting from a local Gibbs state associated to the 
steady state density profile p. In this case, the solution of the hydrodynamic 
equation (|1.2|l is constant in time and equal to p. In particular, the density 
fluctuation field converges, as L ^ oo, to the Ornstein-Uhlenbeck pro- 
cess l|2.2fl with p(-) in place of p^{s,-). Let £ be the differential operator 
defined by 

C = A + X{1 - 2p{x))V (2.3) 

with Dirichlet conditions at the boundary. An elementary computation 
shows that for any smooth function G with compact support on (0, 1) 



6(G) = eo(e*^G) + ^ ^ (Ve^^-'^^^C) (x) V2x{p{x)) W{s,x)dsdx . 



(2.4) 

Since we imposed Dirichlet boundary conditions for £, e'^^G converges to 
as t ^ cx). The second term of the right hand side of (|2.4|) is a Gaussian 
variable with zero mean and variance given by 



t .1 
ds / dx 



(e(*-^)^G) {x)]\{p{x)) 



Integrating by parts the previous expression with respect to x, recalling 
the explicit formula H2.3|) for the operator C and equation (|1.5|l for the 
stationary density profile, we obtain that the previous integral is equal to 

' G{x)\{p{x)) dx - (\e'^G){xfx{p{x)) dx 
Jo 

+ '2J ds dx (e(*~")'^G(x))^ p'{x) . 

Since e'^^G ^ as t oo, the second term of the previous equation 
vanishes as t — > cx3. On the other hand, since p is solution of the stationary 
equation H1.5|l . the operator C is symmetric in L"^ [p' (x)dx) . In particular, 
the last term of the previous equation can be rewritten as 

2 J f ds r dxG{x)[e^^'-'^^G]{x)p'{x) 







which converges, as t — > cxd, to 



J / G{x)[{-C)-^G]{x) p'{x) dx 
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In conclusion, the generalized Ornstein-Uhlenbeck process £,t converges, 
as i ^ cxD, to a Gaussian field with covariance given by 

{eHG),eHF)) (2.5) 

G{x)F{x)x{p{x)) dx + J G{x)[{-C)-^F]{x) p'{x) dx . 



This Gaussian field is also the stationary state for the Ornstein-Uhlenbeck 
process (|2.2f) and the limit of the density fluctuation field under the sta- 
tionary measure. 

An elementary computation permits to rewrite the covariance as follows 



(r*(G),e^*(F)) = (2.6) 
G{x)F{x)x{p{x))dx + f dx f dyF{x)f{x,v)G{y) , 







where / is given by (|1.8|) . Indeed, for F — G, the second term on the right 
hand side of H2.5|) can be written as 

Jsup|2 / dxG{x)H{x)p'{x) - I dxH{x){~£H){x) p'{x)] , 

where the supremum is carried over all smooth functions H vanishing at 
the boundary. Since p is the solution of the equation H1.5|) . an integration 
by parts gives that 

dxH{x){~CH){x)p'{x) ^ [ dx[H'{x)fp'{x) . 



On the other hand, since H vanishes at the boundary, an integration by 
parts permits to rewrite the linear term as 

dxG{x)H{x)p'{x) = -2 dx(^J^ dy p' {y)G{y) - A^H' [x) 

for any constant A. In view of the previous two expressions, it is not 
difficult to show that the supremum in H is equal to 

dx^( r dy p'{y)G{y) - A) ' (2.7) 
with A given by 

la d^ '^Hx) 
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It remains to develop the square and to recall the factor J in front of the 
variational formula to obtain the second term on the right hand side of 
ijT^ with / given by 

Notice that in the symmetric case (when q = \, i.e. A = 0), £ is the 
usual Laplacian A, p' (x) ^ Pb — Pa = — J , and (|2.5|1 becomes 

(f*(G),r*(^)) = C G{x)F{x)x{p{x)) dx-{pb-paf r G{x)[{-A)-'F](x) dx 
Jo Jo 

in accordance with the covariance formula obtained by Spohn in [2] . 



3 Static Approach 

The idea of the derivation is similar to the one in the totally asymmetric 
case In |11| . under the assumption that 

JiPa - Pb) > , 

(which holds, for example, when pa < Pb and g > 1), it was showed that 
the probability dP{{p{x)}) of observing a given macroscopic profile p{x) in 
the steady state can be written in the large L limit as 

dP{{p}) ^ V[{p}] [ e-^^^^P^'^y^^V[{y}] (3.1) 
J{y} 

where the sum is over all positive continuous functions {y{x),0 < x < 1}. 
We will call such a function {y{x)} a path. The expression H3.1|l is written 
there as a path integral, but it is nothing more than the large L limit of a 
sum over discrete paths y{i/L) — and the measure H3.1|) can be simply 
thought as the weight of these discrete paths. 

The function Q{{p{x)}, {y{x)}) was given in JJ, equation (3.22): 

g{{p{x)}; {y{x)}) = -K,{p,, pb) + 2/(0) log + y{l) log ^ ' 



1 - Pa Pb 



dx 







1 - e-^y 



log h plogp + (1 - p) log (1 - p) 

+ (1 - p + y') log (1 - p + y') + (p - y ) log (p - y) 

where K\{pa, Pb) is a normalization constant. The quantity e~^^(l^''l'l^*'l)2?[{j/}]r'[{p}] 
can be thought as the joint probability of the profile {p} and the path {y}. 
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One can rewrite as 

d¥{{p})=V[{p}] f r{{p}\{y})dQi{y}) 

where dQ{{y}) is the probability measure of the positive walk {y} 

dQi{y}) - V[{y}] f e-^e^^J-^^J'^Kp}] ^ 
with Q({j/})= inf g(M,M), i.e. 

Q{{y}) = -K^ipa, Pb) - log4 + y{0) log + y{l) log ^ ~ 



dx 



-log 



1 - Pa 

1 - e-^y 



A 

+ (1 + y') log (1 + y'} + (1 - y') log (1 - y') 



(3.2) 



(as for a given y, the optimal profile is py = ^^^) and ?'({p}|{2/}) is the 
conditional probability of the profile p given the path {y} given by the 
Radon-Nikodym derivative 

r{{p}\{y}) - ^^'^^ 

so 

log(r({p}|M)) ^ 



log4+ / da;{-plogp- (1 -p)log(l -p) 
Jo 



L 

+ (1 + y') log (1 + y') + (1 - y') log (1 ~ y') 
-{1-P + y') log {l~ p + y')-{p- y') log (p - y')} 



The fluctuations of p[x) in H3.1|l have thus two contributions: one com- 
ing from the choice of the path ?/(x), and the other one from the randomness 
of the profile p{x) once the path y[x) is chosen. We shall see that for small 
fluctuations, these two contributions are uncorrelated, leading to Hl.ll|l . 

The optimal path yopt{x), which maximizes the expression ()3.2|) of 
dQ({t/}), is solution of 

2/;pt(0) = 2p,-l 
2/opt(l) = 2/3b - 1 

— ^ = — ^^^^ • (3.3) 
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By expanding (|3.2|l up to the second order, we get for a path y close to 

Vopt- 

y{x) = yopt{x) + 5y{x) 
that its probabihty measure is given by 



exp —L 



dx 



l\2 



2 (i_e-W)2 l-y'^^. 



As yopt{x) > and 6y ^ the condition that 

y{x) = Voptix) + Sy{x) > 
is automatically satisfied. 

The optimal density profile Py{x) for a given {?/} (i.e. the one which 
maximizes ?'({p}|{j/})) is given by 



Py 



l + y' 



(3.4) 



Given the fluctuation Sy{x) = y{x)—yopt{x) of the walk {y}, the probability 
of a small fluctuation of density 5b{x) = p{x) — Py{x) around py is obtained 
by expanding ''({p}|{2/}) up to the second order in Sb and 5y: 



ri{p}\{y})^e^p(-L / dxiSb) 



pt 



(3.5) 



As ?'({/3}|{2/}) does not depend of 6y at order ((5&)^, the choice of (56 — p—py 
is independent of the choice of the fluctuation of the path 6y. 
The total density fluctuation dp is then given by 



5p = 6b + Py — p 
2 



Sb 



where we used (|3.4ll . Thus, by rewriting in (|3.7|) 6y — and 



B 



Vp(i-p) 



/2L 



, and using 



, (EiH), and the fact that 

_ 1 + yopt 

P= TT^ , 



(3.6) 
(3.7) 

Sb = 



(3.8) 



one gets (|1.11() . 
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We now make the link between the small fluctuations and the large 
deviation functional of the density. From (|3.7I) , we see that the probability 
of a fluctuation 5p of order around the optimal proflle p is given by 
the sum over all path fluctuations 5y of the probability of having a path 
fluctuation 6y and a density fluctuation Sb = 6p — ^ around py, i.e. 



dF{{5p}) ^ V[{5p}] / V[{5y}]e 



TT ■, — TTT- 



opt )2 



(3.9) 



In fact, this expression follows directly from (|3.1|l . Integrating it over 5y 
leads to a gaussian form for the density fluctuations 5p of order which 
is another way of writing H3.7() . On the other hand, considering H3.9|l for 
5p small but of order 1 (i.e. of order in L), and performing a saddle 
point evaluation over 5y leads to leading order in 5p to a deviation func- 
tional quadratic in 5p which is identical to the gaussian. The fact that 
this quadratic form of the large deviation functional (for 5p small but of 
order 1 in L) is equivalent to the expression of the gaussian fluctuations 
(for 5p ~ ^) shows that for the WASEP the fluctuations of order ^ can 
be calculated by expanding the large deviation functional to leading order 
around the most likely profile. Mathematically, this is simply due to the 
fact that the saddle point calculation is exact when one deals with gaussian 
variables (here the 5y). 



4 Derivation of (11.151) 

Let us define 



-p'{x)J 

aix) — 



2x{-p{x)Y 



Writing Y{x) ~ ^(0) + f^Y'{s)ds and performing the Gaussian integral 
over Y{Q) in H1.13|l . one obtains the following expression for T{u,v) in 

(4.1) 

2x(p(a;)) Jaa{z)dz 

We show now that when one writes T~^{x, y) as in (|1.15|) . one gets expres- 
sion (|1.8|) for f{x,y). 

Firstly, the symmetry of T^^{x, y) implies the symmetry of f(x,y): 

f{x,y) = f{y,x) . 
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Then, by definition of tlie inverse of an operator, we have 



5{u-v)= I T{u,t)T-\t,v)dt 



(4.2) 



Inserting l|4.1|l and (|1.15(l in (|4.2(l gives the following integral equation for 



f{x,y) _^ 1 



^X{p{x)) a{t)dt 



+ dt a{z)dz I a{z')dz'f{t,y) + 
Jo Jo Jx 

, ©(a; - y)x{p{y)) 



dt j a{z)dz J a{z )dzf{t,y) 

Q{y - x)x{p{y)) 



a{z)dz / a{z')dz' 



a{z)dz I a{z')dz' 



H4.3|) implies that for < ?/ < 1, one has 

/(0,2/) = /(l,y) = 



(4.3) 



(4.4) 



Furthermore, when one applies the operator 9a;{^jfy-} to (|4.3|) . one gets 
the following differential equation: 



df-A^^-fi.,y)-5i.-y)^=, 



(4.5) 



a{x) 

This equation can be written as 

dlf[x, y) + {{2p{x) - l)/(x, y)} + 5{x - y)Jp'[x) = (4.6) 

The 5{x — y) function simply means that the dxf{x,y) is discontinuous in 
x = y, i.e. 



dif{x ,x) - dif{x+ ,x) = Jp'{x) 



(4.7) 



This differential equation with boundary condition (|4.4I) is equivalent to 
equation H1.10|) and its solution is given by H1.8|l (using equation (|1.5|) and 
the symmetry of f{x,y) ) . 



5 Conclusion 

We have proved that the density fluctuations in the stationary state of the 
weakly asymmetric exclusion process with open boundaries are distributed 
like a Gaussian field. We have obtained a simple expression 1)1. 6|l . H1.8|l of 
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the two point function of the weakly asymmetric exclusion process which 
extends the result of Spohn [2]. This correlation function is long ranged 
and is a signature of the non locality of the large deviation function. As for 
the symmetric case, expanding the large deviation functional of the density 
around the optimal profile leads to the right expression. Our results have 
also some similarity with those of the totally asymmetric case (TASEP) 
There too, the fluctuations of density can be written as a sum of two 
terms. However for the TASEP one of the two processes (the Brownian 
excursion) is non Gaussian. It is to be noted that this Brownian excursion 
is not the large A limit of the process Y{x) and, more generally, 

the correlation function does not converge in the large A limit to the 
correlation function of the totally asymmetric case, meaning that the large 
L limit and the large A limit can't be inverted. 
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